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With the relativistic Coulomb wave function boundary condition, the energies, widths and wave 
functions of the single proton resonant orbitals for !" Ne are studied by the analytical continuation 
of the coupling constant (ACCC) approach within the framework of the relativistic mean field 
(RMF) theory. Pairing correlations and contributions from the single-particle resonant orbitals 
in the continuum are taken into consideration by the resonant Bardeen-Cooper-Schrieffer (BCS) 
approach, in which constant pairing strength is used. It can be seen that the fully self-consistent 
calculations with NL3 and NLSH effective interactions mostly agree with the latest experimental 
measurements, such as binding energies, matter radii, charge radii and densities. The energy of 
728;/5 orbital is slightly higher than that of 71d5;? orbital, and the occupation probability of 
the (72s, jo)” orbital is about 20%, which are in accordance with the shell model calculation and 


three-body model estimation. 
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I. INTRODUCTION 


The study of unbound or loosely bound exotic nuclei casts a new light on the discovery 
of nuclear halo structure [1]. The valence nucleons can be readily scattered into single- 
particle resonant orbitals in the continuum. Due to the small binding energy and small or 
no centrifugal barrier conditions, the valence nucleons tunnel out of the potential barrier 
to long distances with an extended density tail to form the so-called halos. Therefore, the 
resonant orbitals in the continuum and the coupling between bound states and the continuum 
near the threshold play an important role in the description of halo phenomena [2, 3]. In 
most calculations, the continuum is replaced by a set of positive energy orbitals without the 
contributions of the widths, determined by solving the Hartree-Fock-Bogoliubov (HFB) or 
relativistic-Hartree-Bogoliubov (RHB) equations in coordinate space and with box boundary 
conditions [4, 5]. Taking into account of the contributions from bound states and resonant 
orbitals in the continuum, the resonant Bardeen-Cooper-Schrieffer approach based on the 
relativistic mean field theory (RMF-rBCS) [6-10] is applied to describe the properties of 
neutron-rich nuclei. 

Recently, neutron halos have been prominently observed in experiments |1, 11-14] and 
widely explored theoretically [15-19]. However, studies on proton halos are still scarce 
because of the existence of the Coulomb barrier [20-22], and even more so for two-proton 
halos [23]. As one of the promising two-proton halo candidates, Ne recently attracts 
experimental [24-26] and theoretical [27-32] attention. The most precise mass and charge 
radii measurements of !’~??Ne have been performed with Penning trap mass spectrometry 
and collinear laser spectroscopy [25] Furthermore, the reaction cross sections for Ne are 
measured by the transmission method and the corresponding density distribution is deduced 
from a modified Glauber-type calculation [26]. 

The improved experimental measurements for !"Ne make it possible to better compare 
with theoretical calculations. Using a Woods-Saxon potential together with Coulomb inter- 
action in three-body model, the ground state properties of "Ne are analyzed [28] in which 
the pairing correlation parameters are adjusted to reproduce the threshold of two-proton 
emission. The root mean square (rms) radius for !"Ne is estimated to be 2.62 fm and the 
s-orbit occupation probability is about 1596. It is also useful to make similar calculations 


for Ne in our fully self-consistent model, including proton resonant orbitals contributions, 


pairing correlations and Pauli block. 
The paper is organized as follows. Numerical details, results and discussion are given in 


Section II. The conclusions and outlook for future work are outlined in Section III. 


Il. NUMERICAL DETAILS, RESULTS AND DISCUSSION 


The single-particle (s.p.) energies Eg and wave functions for bound states are extracted 
from the RMF approach. With the relativistic Coulomb wave functions boundary condi- 
tion [33], the s.p. energies Eg, widths I, and wave functions for proton resonant orbitals are 
extracted from the analytical continuation of the coupling constant (ACCC) [34] approach 
within the framework of the RMF theory. Using Eg, Eg and I, we solve the BCS gap 
equations to obtain the Fermi energy, pairing gap, the occupation probabilities of bound 
and resonant orbitals [10]. With contributions from the relevant s.p. orbitals, the nuclear 
densities are recalculated solving the coupled Dirac equations, as well as equations of meson 
field and photon field. The total binding energy is recovered from the summation of s.p. 
energies, meson energies, pairing energies and center-mass corrections. The next iteration 
starts with new Eg, Eg, I and wave functions for s.p. orbitals extracted from the new 
potential. Therefore, the RMF+ACCC+BCS calculations are numerically self-consistent by 
an iterated way until convergency. The binding energies are precise to 1073 MeV and the 
densities to 1074 fm~*. Details for the neutron case can be found in Ref. [10]. 

The pairing window is opened from all the bound states below the Fermi surface to 
one harmonic oscillator shell above the Fermi surface in our model, including the bound 
proton states 18;/5, 1p3/2, 1p1/2, and single proton resonant states 1d5/2, 28;/2. The constant 
pairing strength G = C/A in the BCS approximation, is chosen by fitting the odd-even mass 


difference extracted from three-point formula for neutron pair, 


A, = s{B(Z.N -1) — 2B(Z, N) + B(Z, N +1)], (1) 


in which binding energies are taken from Ref. [35] except for Ne by latest mass measure- 
ments [25]. For "Ne, C, = 25.5 (26.0) can reproduce A, = 1.8 MeV for NL3 (NLSH) 
effective interaction. By fitting the odd-even mass difference A, — 0.9 MeV extracted from 


three-point formula for proton pair, 


A ECIBUZ d NT EU NEBULA (2) 


C, = 10.5 (11.0) for NL3 (NLSH) effective interaction. Highly excited resonant orbitals with 
large widths, such as 1d3/2 and 1f7/?, are ignored in our calculations because of the minor 
occupation probability. For odd neutrons, the blocking effect is considered in the pairing 
correlations [36]. Note that RMF+ACCC+BCS approach is different from the conventional 
BCS approach, taking only the real resonant orbitals into consideration and kicking the spu- 
rious continuous orbitals off from the continuum [9, 10]. Using this approach, the properties 
of s.p. bound orbitals and resonant orbitals in stable and unstable Ni, Zr and Sn isotopes 
have been well described [10, 37, 38}. 

The theoretical binding energies for "Ne and !°O are calculated by the RMF, RMF+BCS 
and RMF+ACCC+BCS approaches with NL3 and NLSH effective interactions [39, 40]. As 
is well known, the results of the RMF+BCS approach (i.e., the conventional BCS approach) 
change with the box size because of the inclusion of the spurious orbitals. As opposed to 
the RMF--BCS approach, the results of the RMF+ACCC+BCS approach are invariable for 
the specified pairing window, which is from all the bound states below the Fermi surface 
to one harmonic oscillator shell resonant orbitals above the Fermi surface. Therefore, only 
the binding energies obtained from the RMF and the RMF+ACCC+BCS approaches are 
shown in Tab. I. The available experimental binding energy for !"Ne is 112.9276 MeV and 
that for 15O is 111.955 MeV [35]. An asterisk denotes the result of the RMF+ACCC+BCS 
calculation without resonant width contribution. In this case, the resonant orbitals are 
simplified by the discrete positive orbitals without the widths. It can be seen from Tab. I that 
the results without width effect are larger than those with widths contributions because of 
the over evaluation of the pairing correlation, which is also mentioned in Ref. [9]. The binding 
energies in the RMF+ACCC+BCS calculations are 1% larger than experimental data. The 
two-proton separation energy S», for "Ne is 0.93 MeV, and the corresponding theoretical 
result is 0.2 (0.14) MeV for the NL3 (NLSH) effective interaction. The values of S5, in our 
calculations are smaller than data, but those at least tell us !Ne is a weakly bound nucleus, 
which can not be seen from the RMF calculations. The binding energy for !?^O is larger than 
that for !"Ne in the RMF calculations, which means !"Ne is unstable against proton emission 
if pairing correlations are not taken into account. In the RMF+ACCC+BCS calculations, 
the binding energy for PO is less than that for Ne, which means the valence protons are 
weakly bound in !'"Ne because of the attractive pairing correlations. As a consequence, 


pairing correlations are crucial in the existence of the nuclide and need to be taken into 


consideration. The contributions from the resonant orbitals in the continuum to the pairing 
are as important as the bound states near the Fermi surface. The width of the resonant 
orbital has the effect of weakening the pairing gap and decreasing the binding energy. All 
these factors are self-consistently included in the RMF+ACCC+BCS calculations so that 
the reasonable description of the binding energy can be given. 

To show the contributions from the s.p. energies to the nuclear binding energies, the 
single proton orbitals of !"Ne are displayed in Tab. IT, which are calculated by the fully self- 
consistent RMF+ACCC+BCS approach compared with those by RMF approach, for the 
NL3 and NLSH effective interactions, respectively. Note that the positive energies of the s.p. 
orbitals given by the RMF approach are improper because of the box boundary condition. 
The results of the RMF approach are also shown for comparison with the results of the 
RMF+ACCC+BCS approach. It can be seen that the energies of the bound states in the 
RMF calculations are slightly different from those in the RMF+ACCC+BCS calculations 
because the self-consistent potential of the latter is not the same as the former. For protons, 
the bound orbitals 1851/2, 1p3/2, 1pi/2, as well as the resonant orbitals 1d5/? and 2s;/2 with 
narrow widths, have a similar structure for the NL3 and NLSH effective interactions. For 
clarity, the s.p. energies of proton orbitals in '’Ne are plotted in Fig. 1 in potential well and 
barrier, for NL3 effective interaction. The proton potential with the Coulomb barrier V, (red 
curve) corresponds to the potential for s partial wave; V, +1(1+1)/r? (blue curve) refers to 
the potential with the Coulomb barrier and centrifugal barrier for d partial wave. It can be 
clearly seen that the 2s;/? orbit (red line) lies below the Coulomb barrier (J = 0), and the 
1ds;5 orbit (blue line) lies below the Coulomb and centrifugal barrier (l = 2). Since s and d 
orbitals with positive energies lie just above the threshold (E=0) and below the barrier, the 
two valence protons have the probabilities to tunnel the barrier, which are proportional to 
the widths of the resonant orbitals. 

Furthermore, the density distributions for the core nuclide !^O and two-proton nuclide 
"Ne in our model are presented in Fig. 2 in logarithm scale, for the NL3 (solid line) and the 
NLSH (dashed line) effective interactions, respectively. The error region for the experimental 
results is displayed by the gray area [26]. The radii of ^O and "Ne are also shown in Fig. 2. 
It can been seen that the densities from the RMF+ACCC+BCS approach agree well with 
the experimental fit results. The available experimental rms radius of !"Ne is 2.75(7) fm 


and experimental charge radius is 3.042(21) fm. In the RMF+ACCC+BCS calculations, 


the rms radius of "Ne is about 2.86 fm for the NL3 effective interaction, and 2.79 fm for the 
NLSH effective interaction; the charge radius of !"Ne is about 3.15 fm for the NL3 effective 
interaction, and 3.07 fm for the NLSH effective interaction. The relative errors are less than 
5%. 

To investigate detailed contributions to the radii, we plot the occupation probabil- 
ities of the single proton orbitals 1s;/5, lps/5, Ipiy2, 1ds/2 and 2s;/2 in INe by the 
RMF+ACCC+BCS approach in Fig. 3. In our calculations, the occupation probability 
of (251/53)? is about 20%, and that of (1d5/2)? is about 40% for the NLSH effective inter- 
action, which is in accordance with the results in shell-model calculations and three-body 
model [28, 41, 42]. In Fig. 4, the real part of the wave functions for the upper component 
of proton resonant orbitals 1d5/2 and 2s;,/5 in "Ne by the RMF--ACCC--BCS approach 
are respectively plotted for the NL3 and NLSH effective interactions. The extension of the 
wave functions for 1d5/2 and 2s,/? can be clearly seen at r > 3 fm, which give the main 
contributions to the density of this range for !"Ne. 

Experimental studies [24-26, 43] have not been able to confirm the halo structure of 
"Ne, and the corresponding theoretical studies are still controversial. The possibility of the 
halo existence for !"Ne is discussed by the occupation probabilities of the resonant orbitals 
28,/j5 and 1ds5/2. With a three-cluster model, Timofeyuk et al. obtain a (2s; /2)°-dominant 
configuration [44]. Some three-body model calculations suggest almost equal occupation 
probabilities of the (2s1/2)? and (1ds/2)? orbitals [27, 45, 46], but some calculations suggest 
a (1d5/2)? dominance [28]. From the Coulomb mass shift, Nakamura et al. suggest a (2s1/2)? 
dominance [47], but Fortune et al. oppositely suggest a (1ds5/2)? dominance [41]. By cal- 
culating the interaction cross section o, with a Hartree-Fock type wave function and the 
Glauber model [48], Kitagawa et al. also proposed a (1d5/2)? dominance. 

In our present scheme, the occupation probability of (1d5/2)? is larger than that of (25;/5)?, 
because the 72s 1/2 orbital lies slightly above v1d5;? orbital. The inclusion of tensor force 
might alter the present results, which has been discussed for proton 1d5/2—1d3/2 splitting 
decreasing for Ca isotopes in Ref. [50] within the framework of the RMF model, also ap- 
pearing in Ref. [51] within the framework of HFB+SLy5+tensor force model. Moreover, the 
consideration of pairing correlations is important since the filling near the Fermi surface is 
embedded in the continuum so that the density of the energy levels becomes large. In this 


sense, the state-dependent pairing gap may be necessary. Such an effort on this point is in 


process. 


II. SUMMARY 


The energies, widths and wave functions of the single proton resonant orbitals for !"Ne are 
studied by the RMF-ACCC approach with the relativistic Coulomb wave functions bound- 
ary condition. Pairing correlations and contributions from the s.p. resonant orbitals in the 
continuum are taken into consideration by the resonant BCS approach, in which constant 
pairing strength are used. The fully self-consistently microscopic RMF+ACCC+BCS cal- 
culations reproduce the main experimental properties of two-proton halo candidate !"Ne, 
such as matter radii, charge radii, densities with the NL3 and NLSH effective interactions. 
The occupation probability of the proton resonant orbit (725;/5)? is in accordance with the 
results of the shell model and three-body model. The resonant orbitals in the continuum 
play an important role in the pairing correlations as well as those bound states near Fermi 
surface. Our model considers pairing correlations for proton-rich nuclei by the coupling be- 
tween a few narrow resonant orbitals and bound orbitals, which can be safely extrapolated 
to describe exotic nuclei. Since many nuclei are deformed, it is desirable to include the de- 
formation effect to study the properties of the deformed nuclei. Preliminary attempts have 
been made in a deformed Woods-Saxon potential for coupled Dirac equations, which will be 


explained in detail independently. 
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Tab. I: Binding energies of !"Ne calculated by the RMF and the RMF+ACCC+BCS (this work) 
approach with NL3 and NLSH effective interactions, respectively, and compared with experimen- 
tal binding energy (B(!7Ne)—112.928 MeV, B(1?0)—111.955 MeV) [35]. Two-proton separation 
energies $5, are also shown for the RMF and the RMF+ACCC+BCS approach. All energy values 
are in unit of MeV. Asterisk denotes the result of the RMF+ACCC+BCS calculations without 


resonant widths contributions. 


approach Bnrs Byrs — PEgxp. BNLsH PNisH — Prxp. 
Ne RMF+ACCC+BCS 114.50 1.57 113.68 0.75 
RMF+ACCC+BCS* 114.84 1.91 114.01 1.08 
RMF 112.69 —0.24 111.86 —1.07 
150 RMF+ACCC+BCS 114.30 2.345 113.54 1.585 
RMF 114.29 2.33 113.52 1.56 
Sop RMF+ACCC+BCS 0.20 0.14 
RMF —1.60 —1.66 
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Tab. II: Single proton energy orbitals for !"Ne in the RMF and the RMF+ACCC+BCS calculations 


with NL3 and NLSH effective interactions, respectively. All energy values are in unit of MeV. 


vnl; 


281/2 
idis 
1p1/2 
1p3/2 


18172 


NL3 NLSH 


RMF RMF+ACCC+BCS RMF RMF+ACCC+BCS 


E E Tr E E Tr 
0.81 0.73 0.02 1.24 1.17 0.18 
0.49 0.53 0.003 0.61 0.65 0.0026 

—9.12 —8.94 0 —9.23 —9.09 0 
—15.53 —15.36 0 —15.98 —15.84 0 
—35.71 —35.43 0 —36.31 —36.08 0 
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Fig. 1: (Color online) Single proton orbitals in !"Ne calculated by the RMF--ACCC--BCS approach 


for the NL3 effective interaction. 
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Fig. 2: Density distributions in logarithm scale for the core nuclide !°O and !"Ne in the 
RMF+ACCC+BCS calculations with the NL3 (solid line) and the NLSH (dashed line) effective 
interactions, respectively. The gray region corresponds to the error region for experimental re- 


sults [26]. 
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Fig. 3: (Color online) Occupation probabilities of the proton s.p. orbitals 1s1/2, 1p3/2, 1p1/2; 1d5/2 
and 2s1/2 in 1ITNe by the RMF+ACCC+BCS approach. 
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Fig. 4: (Color online) Upper wave functions for proton resonant orbitals 1ds5/2 and 2s; /2 in "Ne 


by the RMF+ACCC+BCS approach. 
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